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Abstract 



■ We investigate the possible occurrence of a positive cosmic acceleration in a spatially averaged, expand- 

| ing, unbound Lemaitre-Tolman-Bondi cosmology. By studying an approximation in which the contribu- 

tion of three-curvature dominates over the matter density, we construct numerical models which exhibit 
I/*") ' acceleration. 

o 

1 Introduction 

43 ; 

^p H , We live in an inhomogeneous Universe, whose exact and complicated dynamics is described by Einstein's 
equations. It is generally assumed that when the spatial inhomogeneities are averaged over, the resulting 
Universe is described by the standard Friedmann equations for a homogeneous and isotropic cosmology. 
However, as is known pQ, since Einstein equations are non-linear, the averaging over the inhomogeneous 
matter distribution will in general not yield the solution of Einstein equations which is described by the 
Friedmann-Robertson- Walker (FRW) metric. There will be corrections to the FRW solution, which could 
be small or large, and which could in principle lead to observational effects indicating a departure from 
standard FRW cosmology. 

The possibility that the observed cosmic acceleration j^j is caused by the spatial averaging of the 
observed inhomogeneities, rather than by a dark energy, has been investigated and debated in the literature 
[HI ISIS]- A systematic framework has been developed for describing the dynamics of a modified Friedmann 
universe, obtained after spatial averaging [B]. It has been suggested that, within the framework of standard 
cosmology with cold dark matter initial conditions, an explanation of the acceleration in terms of averaged 
inhomogeneities is unlikely to work [7]. However, it is perhaps fair to say that the matter cannot be treated 
as completely closed, and further studies are desirable 8 . 

The Lemaitre-Tolman-Bondi (LTB) cosmology [H], being an exact solution of Einstein equations for 
inhomogeneous dust matter, provides a useful toy model for investigating the possible connection between 
acceleration and averaging of inhomogeneities. Various authors have examined different aspects of the 
model in this regard. The redshift-luminosity distance relation in an LTB model and its possible connection 
with cosmic acceleration, or the lack of it, have been studied by Celerier ^01 5 Ames et al. and by 
Vanderveld et al. [TSJ, Sugiura et al. ^3], Mustapha et al. ^H, Iguchi et al. ^5]. Nambu and Tanimoto 
give examples of cosmic acceleration after averaging in an LTB model consisting of a contracting region 
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and an expanding region. Other works which study cosmic acceleration in LTB models are those by Moffat 
[T7| . Mansouri ^S], Chuang et al. \T3\, Rasanen |2U] and Apostolopoulos et al. [2*T] . 

It has sometimes been suggested in the literature that both an expanding and a contracting region 
are needed for acceleration. In the present paper we will address a question which does not seem to have 
been addressed in the above-mentioned works: can spatial averaging in a universe consisting of a single 
expanding LTB region produce acceleration? We show that the answer is in the affirmative. We do this by 
considering a low density, curvature dominated unbound LTB model in which the contribution of matter 
density is negligible compared to the contribution of the curvature function. Further, we concentrate on 
the late time behaviour of such a model. As a result of this proposed simplification, the calculation of the 
acceleration of the averaged scale factor becomes relatively simpler and conclusions about acceleration can 
be drawn, for specific choices of the energy function. 

In Section [2] of the paper we recall the effective FRW equations, resulting from spatial averaging in a 
dust dominated spacetime. In Section El we discuss spatial averaging for the marginally bound LTB model 
and point out there can be no acceleration in this case. The unbound LTB model is investigated in Section 
|1J in the approximation that the spatial curvature (equivalently, the energy function) dominates over the 
dust matter density, and numerical and analytical examples of acceleration are given. 



2 Averaging in Dust Dominated Spacetime 

For a general spacetime containing irrotational dust, the metric can be written in synchronous and comoving 
gauge 1 , 

ds 2 = -dt 2 + hij(x, t)dx i dx j . (1) 

The expansion tensor ©*• is given by G*- = (l/2)h tk hkj where the dot refers to a derivative with respect 
to time t. The traceless symmetric shear tensor is defined as <r*- = G*- — (Q/3)<5*- where G = G^ is the 
expansion scalar. The Einstein equations can be split [H] into a set of scalar equations and a set of vector 
and traceless tensor equations. The scalar equations are the Hamiltonian constraint (|2aj) and the evolution 
equation for Q (|2b|) . 

(3) 7£ + -Q 2 - 2a 2 = WirGp (2a) 
3 

{3) TZ + G + Q 2 = UirGp (2b) 

where the dot denotes derivative with respect to time t, ^TZ is the Ricci scalar of the 3-dimensional 
hypersurface of constant t and a 2 is the rate of shear defined by a 2 = (I/2)<7*cr?. Eqns. (|2llj) and (|2b|) can 
be combined to give Raychaudhuri's equation 

G + -Q 2 + 2a 2 + 4vrGp = 0. (3) 
3 

The continuity equation p = —@p which gives the evolution of p, is consistent with Eqns. (|2a() . ()2b|) . We 
only consider the scalar equations, since the spatial average of a scalar quantity can be defined in a gauge 
covariant manner within a given foliation of space-time. For the space-time described by (^Q), the spatial 
average of a scalar *S?(t, x) over a comoving domain T> at time t is defined by 

= tt I d 3 xVh^ (4) 
Vv Jv 

where h is the determinant of the 3-metric hij and Vp is the volume of the comoving domain given by 
Vd = f v d 3 xVh. 

1 Latin indices take values 1..3, Greek indices take values 0..3. We set c = 1. 
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Spatial averaging is, by definition, not generally covariant. Thus the choice of foliation is relevant, 
and should be motivated on physical grounds. In the context of cosmology, averaging over freely-falling 
observers is a natural choice, especially when one intends to compare the results with standard FRW 
cosmology. Following the definition Q the following commutation relation then holds [H] 

(f) c -(f)D = Wj)-(f)D(e> B (5) 

which yields for the expansion scalar 

(ey v - (e) v = (e 2 ) v - (0)1 . (6) 

Introducing the dimensionless scale factor av = (Vv/Vvin) 1 ^ normalized by the volume of the domain T> 
at some initial time ij n , we can average the scalar Einstein equations (|2a|) . (|2b|) and the continuity equation 
to obtain [H] 

(6)2, = 3^ , (7a) 
av 

3 (^) 2 - 8nG(p)v + \(R)v = , (7b) 

3^j+4TTG(p)v = Qv, (7c) 

(P)v= -(Q)v( P )v= -3—{p)v. (7d) 

av 

Here {TZ)y, the average of the spatial Ricci scalar ( 3 'TZ, is a domain dependent spatial constant. The 
'backreaction' Qv is given by 

Qv^l((Q 2 )v-(@)l)-2(a 2 )v (8) 



and is also a spatial constant. The last equation (|7d|) simply reflects the fact that the mass contained in a 
comoving domain is constant by construction : the local continuity equation p = — Qp can be solved to give 
p\[~h = po\fho where the subscript refers to some arbitrary reference time to- The mass Mv contained 
in a comoving domain T> is then Mv = f-p pyf~hd?x = j v po\/hod i x = constant. Hence 

My 

\Piv = 77 T ( 9 ) 



which is precisely what is implied by Eqn. (J7dJ). 

This averaging procedure can only be applied for spatial scalars, and hence only a subset of the Einstein 
equations can be smoothed out. As a result it may appear that the outcome of such an approach is severely 
restricted, and essentially incomplete due to the impossibility to analyse the full set of equations. However 
one should note that the cosmological parameters of interest are scalars, and the averaging of the exact 
scalar part of Einstein equations provides the requisite needed information. A more general strategy would 
be to consider the smoothing of tensors, which is beyond the scalar approach that certainly provides useful 
information, albeit not the full information. 

Equations ()7b|) . (|7c|) can be cast in a form which is immediately comparable with the standard FRW 
equations [22]- Namely, 

av AttG „ . , 

-^(Pefr + 3P c fr) (10a) 



av 



a v \ 8?rG 

— = ~^-Peff (10b) 
av J 3 
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with p e ff and P e fj defined as 



/ v Qv (K)v p Qv (n) v 



A necessary condition for (|ll)a|) to integrate to l)10b|) takes the form of the following differential equation 
involving Qx> and (TZ)t> 

Qv + 6—Qv + (n) v + 2^(K) v = (12) 
ax> ax> 

and the criterion to be met in order for the effective scale factor ax> to accelerate, is 

Qv > 4ttG( P )v ■ (13) 



3 The LTB Solution 

The system of equations (jlUaj) . ()10b|) and lfT2*|) is only consistent, it does not close. For a completely general 
spacetime with dust, therefore, it is not possible to proceed with the analysis without making certain 
assumptions about the form of the functions Qx> and (TZ)x> |H1E]- For this reason, it becomes convenient 
to work with the LTB metric, an exact solution of the Einstein equations which is a toy model consisting 
of a spherically symmetric inhomogeneous dust dominated spacetime. In this section, we describe the LTB 
solution and apply to it the averaging procedure described above for the simplest, marginally bound case. 
In the next section we extend the analysis to the unbound LTB solution. The LTB metric for pressureless 
dust is given in the synchronous and comoving gauge, by 

^2 = _ dt 2 + R' 2 (r t) 2 + r2 ^ t) ^ 2 + s . n2 M ^ (M) 
1 + 2E[r) 



The Einstein equations simplify to 



^-w$= Eirh (15a) 

gfrlwM) - (15b) 

Surfaces of constant r are 2— spheres having area AirR 2 {r,t). p(r,t) is the energy density of dust, while 
E(r) and M(r) are arbitrary functions that arise on integrating the dynamical equations. Solutions can 
be found for three cases E(r) > 0, E(r) = and E(r) < 0. We will restrict our attention to models in 
which E(r) has the same sign, for all r. The solution for E(r) = (the marginally bound case) has the 
particularly simple form 

1 /3 

R(r,t)=(^^-j (t-t (r)f\ for£(r) = 0. (16) 

Here to(r) is another arbitrary function arising from integration. The solution describes an expanding 
region, with the initial time ij n chosen such that t > ti n > to(r) for all r. For the other two cases, the 
solutions can be written in parametric form 

R= G 2 £^ (coshf?-l) ; t-io(r) = (sinh 17-17), < t? < 00 , fbr£(r)>0. (17a) 
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i?= ^1^(1 -cost/) ; t-t (r) = — ^^—(rj-smr,) , < 77 < 2tt, for£;(r)<0. (17b) 

In the unbound case (E(r) > 0), R(r,t) increases monotonically with t, for every shell with label r. In the 
bound case (E(r) < 0), R(r,t) increases to a maximum value R m ax{f) for each shell r and then decreases 
back to in a finite time. 

In all cases, there are two physically different free functions, although three arbitrary functions E, M 
and to appear. One of the three represents the freedom to rescale the coordinate r. We use this freedom to 
set R(r,ti n ) = r. To completely specify the solution, we specify the initial density pi n (r) and the function 
E(r). This specifies M(r) = Air J Q r pi n (r)r 2 dr (which in the marginally bound case is interpreted as the 
mass contained in a comoving shell), and to(r) can be solved for using Eqns. JIHJ), (fTTajl or Ip7h"j) as the 
case may be, at time t = ti n . 

3.1 Averaging the LTB Solution 

The quantities defined in Sec. 12 can be computed for the LTB metric of Eqn. (|14|) . The averages are 
computed over a spherical domain of radius r%>, centered on the observer. Other choices of the averaging 
domain will possibly yield different results, however, the choice of a spherical domain seems natural for the 
spherically symmetric metric of Eqn. (|14|). For clarity, we suppress the r and t dependences of the various 
functions in the following 

, n> R'R 2 , 4vr f 13 R 2 R' + 2RRR' , V v „, f rv M 1 , . . 

V v = 4tt - dr ; O p = — / — dr = -f- ; M v = - dr 18 

Jo VT+2E x ' Vv h V1 + 2E Vv Jo VT+2E v ' 

where a prime denotes a derivative with respect to r. Note that Mp 7^ M{rv) if E 7^ 0. Only in the 
marginally bound case is the function M(r) identified with the mass contained in the shell with label r. It 
is convenient to work with the combination (2/3)(0 2 )x> — 2(a 2 }v rather than evaluate the average rate of 
shear (cr 2 )v separately. We define this to be Cv and obtain 

2 2 2 8vr r 13 2RRR' + R 2 R' 



Cvs - mv -2 { ^ D = -J g ir. (19) 



3.2 The marginally bound case - vanishing backreaction 

The results of the previous subsection hold for all classes of the LTB solution, provided the averaging 
domain is spherically symmetric about the center. Now consider the marginally bound case E(r) = for 
all r. The algebra in this case becomes very simple, and the backreaction can be computed analytically. 
We will show next that for a single domain with E{r) = throughout, the backreaction Qv is, in fact, zero. 
Also, the average spatial curvature (TZ)v is zero (which is expected by inspection of the metric (fTljl if we 
note that a spatially uniform initial density profile in the LTB solution in this case yields the corresponding 
FRW solution). As described earlier we have 



r 3/2 / 2 

M(r) = 4vr / Pin (r)r 2 dr ; t (r) = t in - — J . (20) 







GM(r) 



To obtain the second equation we have used Eqn. (|16|) at time t = t% n with the condition R(r,ti n ) = r. 
Some algebra then yields the following results 



^(«j-'- (( _ (o(r)) -v 3 (21a) 



1/3 



5 



2 / 9GM(r) 

3 V 2 

■ l= 2 ( 9GM(r) 
9 1 2 



1/3 



1/3 



(t-t (r)) 



-1/3 



-4/3 



V© = 6vrG'M I ,(i-toM) z 



M'(r) 
2M(r] 

'M'(r) 
M(r) 



ax»(i) 



(t-«d(r))-t{,(r) 



(t-t (r)) + t' (r) 



tojry) \ 



2/3 



to(rv)J 



(21b) 



(21c) 



(21d) 



Cv 



(21e) 



-to(r©) ' 3(t-toM)" 

Since the solution was constructed assuming t > io( r ) f° r an r J Eqn. (|21cl|) immediately shows that 



dx> < and hence acceleration is not possible in this case. Further, Eqn. (|21e|) shows that Qx> = Cx> — 
(2/3) (0)?) = 0. Thus the backreaction term vanishes for a region described by the marginally bound LTB 
solution. This result is not unexpected. We note that mathematically, the General Relativistic equations 
(|15a|) and (|15b|) describing the evolution of a spherical dust cloud are identical to the corresponding 
Newtonian equations. It has been shown by Buchert, et. al. |23| that the backreaction Qx> in a spherically 
symmetric Newtonian model of dust, must vanish. Further, we note that in the marginally bound case, 
the mathematical expressions for the averaged quantities defined earlier coincide with their corresponding 
Newtonian analogues. Hence, for the fully relativistic (marginally bound) case also, the backreaction must 
vanish. 

The spatial Ricci scalar and its spatial average for a general E(r) are given by 
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E + 



E'R 



16vr 



r-p 



d 
dr 



(ER) 



V1 + 2E 



dr 



(22) 



which shows that ^'IZ and hence {H)t> vanish in the marginally bound case. This is consistent with the 
requirement of Eqn. (|12j) . 



4 The unbound LTB solution 

Since the solution with zero spatial curvature fails to produce a non-trivial backreaction, we consider next 
the opposite extreme - a curvature dominated solution in which the contribution to the Einstein equations 
due to matter is much smaller than that due to spatial curvature. Before describing the construction of 
such a solution, we present a general treatment of regularity conditions which an unbound LTB model 
must satisfy. 



4.1 Regularity conditions on unbound LTB models 

Consider the class of unbound LTB models given by (|17a|) . The functions M(r) and E(r) are to be specified 
by initial conditions at t = ti n , and the choice of scaling R(r,ti n ) = r fixes to(r) as 

GM(r) . , . .. , 2E(r)r . . 

to{r) = U n - 3/2 (sinh7fa(r) - fft w (r)) ; cosh r) in (r) - 1 = QM ^ ■ (23) 

The regularity conditions imposed on this model, and their consequences, are as follows 
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• No evolution at the symmetry centre: 

This is required in order to maintain spherical symmetry about the same point at all times, and 
translates as R(0, t) = for all t. The right hand side of Eqn. (|15aj) must therefore vanish in the 
limit r — ► 0. Since the functions involved are non- negative, we assume that we can write 

E(r -» 0) ~ r 5 , 6 > ; M(r -► 0) ~ r a ; R(r -»■ 0, i) ~ r^/(t) , a > /? > . (24) 

Consistency requires /? to be constant, and our scaling choice further requires (3 = 1. We do not 
require the exponents 5 and a to necessarily be integers. 



• No shell-crossing singularities: 

Physically, we demand that an outer shell (labelled by a larger value of r) have a larger area radius 
R than an inner shell, at any time t. Unphysical shell-crossing singularities arise when this condition 
is not met. Mathematically, this requires 

R'(r, t) > for all r, for all t. (25) 

• Regularity of energy density: 

We demand that the energy density p(r, t) remain finite and strictly positive for all values of r and t. 
Combining this with Eqns. ()15b|) and (|25|) gives (assuming that R' is finite for all r and since (3 = 1) 



lim p(r, t) = finite a -I- 2/3 = 0^ a = 3 . (26) 

) — >o 

• iVo trapped shells: 

In order for an expanding shell to not be trapped initially, it must satisfy the condition r > 2GM(r). 
Near the regular center, this condition is automatically satisfied independent of the exact form of 
M(r), since there M ~ r 3 . 



Consider now the function to(r) given by (|23jl . By observing the behaviour of the functions (coshr/j n — 1) 
and (sinhr^ — rji n ) for values of 5 equal to, less than, and greater than 2, it is easy to check that to(r) is 
finite at r = for all values of 5. However, this involves the assumption that M(r) is positive for r ^ 0. 
In the limit of M — > for all r, we find 

V2E(t-t (r)) ; t„(r) ^ t m --^=. (27) 

Although now, in the limit r — > 0, to( r ) is finite only when <5 < 2, it will turn out that the integrals involved 
in the averaging procedure are insensitive to the behaviour of to(r) in the r — > limit, and remain well 
defined for all positive values of 5. The expression for ^TZ in (|22[) indicates that the spatial Ricci scalar 
diverges as r — ► unless 5 > 2. However, we note that the spatial Ricci scalar is not a fully covariant 
quantity and depends on our choice of time slicing. The /cmr-dimensional Ricci scalar, obtained after 
taking the trace of the Einstein equations as ^TZ = 87rGp(r,t) is finite at the origin irrespective of the 
behaviour of E(r). It is interesting to see how this cancellation occurs. We have 

On using the Einstein equation (|15a|) we obtain 
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which neatly cancels the contribution from ^1Z, leaving precisely 8irGp(r, t) after applying the second 
Einstein equation (|15bj) . Hence the 4— dimensional Ricci scalar does not impose any further restrictions 
on the form of E(r). The fact that the origin is well behaved can also be seen from the behaviour of the 
Kretschmann scalar, given by 

(A^uvtrn G 2 M' 2 G 2 M M' G 2 M 2 , s 

( )n * ° - 12 w 32 ~mT + 48 ^ ■ (30) 

A condition on the value of 5 is obtained, however, by the regularity of the energy density p(r,t), which 
assumes that R'(r, t) is not only positive, but also finite for all r and t. Equation (|27|) shows that unless 
5 = 2, R' either diverges or vanishes at the center, violating this regularity condition. 

4.2 Late time solution and curvature dominated unbound models 

The function R(r,t) is an increasing function of time in all the unbound models described by (|17aj) . The 
Einstein equation (|15aj) then shows that for sufficiently late times t S> U n , neglecting the term involving 
1/R, all unbound models have the approximate solution given by (|27jl. If on the other hand, we start 
with a model which satisfies GM{r) j \rE{r)) <C 1 for all r, then since our scaling assumes that R = r at 
t = ti n , w e will have GM(r)/(R(r,t)E(r)) <C 1 for all r and for all t > and (|27jl is then an approximate 
solution at all times, the approximation becoming better as t increases. To make this idea more precise, 
consider the closed form expression for t in terms of R obtained by integrating Eqn. (|15aj) 24 

t-^^lS^'^M) ■ n.) = id + ^-^-'(^). (3D 

Hence, imposing R{v, tm) — ^ hcivc 

r 3/2 / Er \ 

* (r) -*-e5£F F (<Hf)- (32) 

Let us write GM(r) / E{r) = eGM (r)/E(r) = eg(r) where e is a dimensionless positive number whose value 
we can control. This relation also defines the functions M(r) and g(r). We can rewrite (J31[) for e <C 1 as 

- W)) = {l + i| in (|) - | (in2 - i) + ((|) 2 ) } . ,33, 

Here 0(x 2 ) represents a power series beginning with a term of order x 2 , and we have used a binomial 
expansion in eg(r)/R and the asymptotic expansion for the inverse hyperbolic sine given by (as x — > 0) 

mi 

sinh" 1 = In 2 - In x + ^x 2 + O (x 4 ) . (34) 

The terms in Eqn. (|33|) involving e vanish as e — ► 0, although the expression in (|33() cannot be inverted 
to get R = R(r,t), due to the presence of the logarithm. We can, however, make the following statement. 
Provided the function g(r)/r is finite for all values of r, then given any starting time ij n , we can choose e 
small enough that the terms involving e on the right hand side of Eqn. (|33|) are all negligible compared 
to unity. Then, since R increases with time, these terms will always be negligible compared to unity. 
Alternatively, given some eg/r = GM/{Er) which is finite for all r, one can always wait for a sufficiently 
long time, and find that the e dependent terms become smaller compared to unity. In this case, we need 
not even assume that e is small. It is in this sense that the approximation involved in writing the equations 
in (|27|) becomes better as t increases (with the caveat that if e is not small, then io( r ) must be given by 
the full expression (|3*2*)) and not the approximation of (127(1 ). This shows that the first of equations is 



8 



the correct late time solution for all unbound models, with the second being a good approximation when 
e is small. The condition that for some e > 0, g(r)/r be finite for all r, and in particular as r — > 0, implies 
that 6 < 2 where 5 is defined in (|24[). This is not inconsistent with the requirement 5 = 2 imposed by the 
criterion of regularity of energy density. 

Consider now a model which begins with negligible matter (e — ► 0) and in which we have waited for a 
sufficiently long time (t 3> ti n ). Eliminating to( r ) from (|27|) the approximate solution becomes 



R(r, t) = v / 2E(r) (t - X t t in ) + X r r (35) 

where we have introduced two placeholders A r and At which will remind us of the relative magnitudes of 
various terms. We will ultimately set A r = Xt = 1. Substituting for R in the expression for the domain 
volume Vd in Eqn. (|18|) . we find 

Vp = (t— XttinflE + X r (t — Xtti n ) 2 lEr + Xf,(t — Xttin)X Er 2 + A 3 Z r 2 (36) 

where we have defined the domain dependent integrals 

A n> V2EE' , A r-° (r-2E)' , 

1 E = 4vr / -. dr ; l Er = 4vr / \ 1 dr 



o VT+2E ' Jo VT+2E 

it, | r 2 • \/2E ) r r v r 2 



x —^l ^mt dr ■ ^= 4 V„ 7rm dr - (37) 

The sum of the exponents of A r and A^ in each term in ()36|) indicates the relative order of that term with 
respect to the leading t 3 term. This approach of treating some terms as small compared to others is valid 
since the various integrals which multiply the powers of t, are all finite and non-zero. 

We note that the solution in Eqn. (|27|) actually corresponds to Minkowski spacetime, since in this 
limit the matter content has been neglected. The corresponding Riemann tensor and Kretschmann scalar 
are exactly zero. The constant time three-spaces are hypersurfaces of negative curvature, with the three- 
curvature being determined by the function E(r). The 'FRW limit of this solution is in fact the Milne 
universe; the solution (|27|) could hence be thought of as the 'Tolman-Bondi' type generalization of the 
Milne universe. For our purpose, it is not a problem that the solution describes Minkowski spacetime - we 
know from the initial conditions that dust matter is present, only its density is negligible compared to the 
curvature term. The form of the solution then allows us to easily determine if the average scale factor ax> 
undergoes acceleration. We demonstrate this with explicit examples in the next subsection. Subsequently, 
we argue that if a small amount of matter is introduced, so as to introduce departure from Minkowski 
spacetime, the sign of the acceleration of a-jj is preserved. 

4.3 Condition for late time acceleration 

The expression for the volume Vx>(t) in (|36|) allows us to determine the late time behaviour of the effective 
scale factor ax>(t) = (Vz>(t)/Vx>(£j n )) x ' 3 . Using a binomial expansion for t 3> ij n in l|3fi|). we get 

3 S-|-5(|) 2 ^(^-3fe^) 2 ) + ° (3) <38) 

where O (3) represents terms involving A™ A", i.e. containing (l/t m+n ) with m + n > 3. 

We see that the generic late time (i.e. t — > oo) behaviour of the unbound models under consideration 
is '&d — > 0, and that deviations from zero are small, being a second order effect. Whether the approach to 
dx> = is via an accelerating or decelerating phase, depends upon the relative magnitudes of the domain 
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Figure 1: The function V defined in the text, for the power law models described by 2E(r) = r s . V is positive 
everywhere except along 6 = 2. 

integrals involved. A sufficient condition for an unbound model with negligible matter to accelerate at late 
times, is 

V =_l Er2 --L(l Er f >0 . (39) 
6±e 

To proceed further we need to specify a particular model. As an explicit example of models admitting 
acceleration, we consider the power law models characterized by 2E(r) = r s , for all r, in some units. (At 
present we are only demonstrating the existence of such models, and shall therefore not worry about the 
physical scales involved.) Keeping in mind the discussion of Sees. 14.11 and 14.21 we must strictly speaking 
only consider the model with 5 = 2. The models with 5 > 2 cannot be considered at all, since they violate 
the conditions assumed in Sec. 14.21 which justified the approximation in Eqn. (|27|) . The models with 
5 < 2 on the other hand, contain a Ricci scalar that diverges and a matter density that vanishes at the 
center. Despite these pathologies, we display the results for the models with 5 < 2 as an existence proof of 
acceleration using this very simple parametrization. Although it is possible to obtain analytical expressions 
for the integrals in (|37[) in terms of the incomplete Beta function, it serves our purpose much better to 
numerically evaluate the integrals for various values of 5 and rx>, and plot the function V defined in (|39|) , 
The results are shown in figQ Note that V vanishes along the line 5 = 2, but is positive everywhere 
else in the region plotted, and that the positivity of dx>/ax> at late times is independent of the size of the 
domain rx>- We have therefore obtained a continuous range of parameter values (5, r©) which admit late 
time acceleration. 

In order to demonstrate that the acceleration obtained above is not an artifact of the singular behaviour 
of those models, we construct another set of models which show late time acceleration, and in which the 
spatial Ricci scalar ^TZ remains finite everywhere. Consider the models characterized by the energy 
function 

2 

2E(r) = — — ; a > (40) 
w l + r a ' y J 

where we have again used arbitrary units. Since a > 0, the r — > behaviour of these models is 2E ~ r 2 , 
which satisfies the regularity conditions of Sec. 14. II and keeps ^7Z finite at the origin. Also, keeping a < 2 
ensures the 'no shell-crossing' criterion of Sec. 14.11 The function V/Ie for these models, which controls 
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Figure 2: The models described by 2E{r) — r 2 /(l + r a ). (a) The scaled function V/Ie- (b) V /Ie plotted against 
a for specific values of rp. 



the magnitude of the late time acceleration, is shown in fig. against r-p and a. For clarity, in the 
second panel we have shown V /Ie against a for specific values of rp. Again, we find that V /Ie is positive 
everywhere in the region shown, and hence the models show late time acceleration for all allowed values of 
a and rp. As an example, we plot the evolution of the dimensionless quantity q-p defined by 

qv = rg— = 2 - 3 — x 2 (41) 



a v 



for various fixed values of a and rp. The results are shown in fig. |31 We have used units in which ti n = 1, 
and have displayed the evolution for times t > 100 ti n . As mentioned earlier, a potentially contentious issue 
is that in all of the calculations above, we have actually set e = 0. Since the function to( r ) approaches 
its approximation in Eqn. (|27|) in a continuous fashion as e — > 0, we expect that models with a non-zero 
but small matter density will also exhibit the same qualitative late time behaviour as the ones above. To 
demonstrate this, we consider the leading corrections to the function to(r) in the presence of a small but 
non-zero e. We are still assuming the late time limit so that the e dependent terms on the right hand side 
of Eqn. (|3.3|) can be neglected (more precisely, we treat both e and g/R as small quantities). First, let us 
rewrite Eqn. ()32j) as 

$8 ^ *)-(^r f ©='{ i+ iK?)^(?)}' <-> 

It is easy to check that in the late time limit, the expression for volume Vx> becomes 

Vp = (t— XtUnflE + K(t — XtUnflEh + Aj!(t — X t ti n )I Eh 2 + A^Z ft 2 (43) 

where Ie is the same as defined in (|37|) ■ and the remaining integrals are defined analogously to those in 
(EH), 

, r° (2EhY , ^ , r c (h 2 V2E)' , t f rv h 2 ti , 
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Figure 3: Evolution of qv in the models with 2E(r) = r 2 /(l + r a ), plotted for (a) three values of a with rv = 250, 
and (b) three values of td with a = 1. 



In arriving at Eqn. (|43|). we have neglected terms involving (g/i?)(e In e), e{g/R)\n{g/R) and terms of 
order 0(eg/R) on the right hand side of Eqn. (|3*H|) . In order to proceed as before, we further assume that 
these leading order corrections are smaller than the terms of order A 2 , coming from the binomial expansion 
of Vx> in Eqn. (|43|) . This is essential in order to be able to make a statement analogous to (|39|) . and can 
be ensured by choosing e small enough, without setting it exactly to zero. The condition for late time 
acceleration in this situation becomes 

V h = l Eh 2-^-{l Eh f>Q. (45) 

For e < e -1 , we have < e < — elne < 1, and the leading order terms in the expansion of h(r) contain 
(elne) (assuming that the r dependent coefficients are well behaved for all r). On expanding the integrals 
in Q45[) to this leading order, we find for the function Vh, 

V h = V + (-elne)P (elne) + O (e, (elne) 2 ) 

n J VTT2E 

where M is defined by M = eM, V is defined in Q39j) and the integrals I^ ne> and give the leading 

order corrections to2"^ r andX S) ,2 respectively. The function Vh of (|45|) replaces T 7 in Eqn. (|38|) . This shows 
that a non-zero e brings in an additional correction to a-D/av which is of order A 2 (elne). We have already 
neglected terms of order (g/i?)(elne), and since (g/R) is small because t is large, we should therefore also 
ignore terms of order A r (elne), At (elne) and higher. Hence the correction given by pi^e) ( provic ied it 
is finite), should be ignored. We therefore see explicitly that within the late time approximation, we can 
always have a non-zero but small enough e which does not affect the sign of the acceleration at the leading 
order. 

A rigorous argument to demonstrate acceleration in models with non-zero matter in general would, of 
course, require a complete numerical evolution of the LTB equations, and this work is in progress. Our 
semi-numerical analysis, however, throws open the possibility of a large class of models which may show 
late time acceleration after averaging. 

Our results for the zero matter limit also point to the need for caution in interpreting acceleration - by 
suitable choices of the energy function E(r) one could obtain acceleration, even though in this limit the 



p(elne) 



^ 7- T 01ne) 



~ 1JJI ~ h,r Er 2 

r v (2GMr/V2E~y 



4-7T 



dr 



(46) 
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spacetime coincides with Minkowski spacetime. A similar demonstration was earlier given by Ishibashi and 
Wald :5j. They showed that by suitably joining two negative curvature slices (hyperboloids) in Minkowski 
spacetime one can construct a spatial region which exhibits acceleration. In contrast though, the slicing we 
choose is physically motivated, so as to coincide with the FRW slicing when matter is included. Our own 
interest of course was in demonstrating, by adding matter beyond the Minkowski limit, that acceleration 
is possible in a single expanding LTB region. 



4.4 An analytical example (r = excluded) 

In this subsection we will follow a slightly different approach and try to construct an accelerating model 
from purely analytic arguments. We begin with a domain in which to(r) > for all r. We now use the 
approximate solution (|27[) in the expression for volume in (|18[) and keep the integration limits unspecified 
as n, V2, i.e. Vp = 47r R'R 2 + 2Edr. At late times t, by treating io( r ) m its entirety as a small 
quantity compared to t, we find 

*> « (47) 

where the constant of proportionality is positive, e dependent terms have been neglected and we have 
defined the integrals 2 

1 r m 3 ' 2 )' dr ; l2 = _ r ((W)' dr ; l3 = r (W)' ir . m 



3j n y/l + 2E J n VI + 2E J ri + 

Note that ax> is still a second order quantity at late times. The integral Zi, which is essentially the same 
as Ze defined in (|37[). is of the form J x 1 / 2 /y/l + xdx. This can be evaluated exactly and X\ is positive 
provided E{r) is an increasing function of r, which we henceforth assume. Let us also assume that to(r) 
is an increasing function of r, this requires E(r) to increase faster than r 2 . Although this is not consistent 
with our arguments of Sec. 14.21 notice that simultaneously requiring to > and t' > places a restriction 
on the minimum value that r can take. For example, if 2E = (r/ro) 171 , m > 2 then these conditions imply 
(using Eqn. (j2*7jl ) that r > r\ > {To l /t 2 n ) 1 '^ m ~ 1 > . Since the origin is necessarily excluded, the function 
GM/(Er) will remain finite for all allowed values of r, even though E(r) rises faster than r 2 . So the 
assumptions on which the later arguments were built, are not violated. 

The condition for positivity of a© at late times is now V = — T 2 /3 > 0. Our assumptions above 
ensure that the integrands of 1%, I2 and X3 are all positive. We denote 

tin ti n (2E)^/ 2 

The positivity condition can be written as 



1 f dr] { I" tUfp2) ' dr] > { f' 2 tmifP) ' drV (50) 
, n 3vT + 2£ j\J ri VT+2E J \J n V3VT+2E J 

A set of sufficient conditions for this relation to hold, is 

f>(fp)' ; (fp 2 )'>(fp)'. (51) 

Some straight-forward algebra reduces these conditions to 

— > and — > (52 

/ 1 ~P J l-p p 



Factors of 4-7T have been absorbed in the proportionality constant. 
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respectively. By definition, < p < 1 since we require to > for all r. This implies (2p — l)/p < 1 for all r, 
so if the first relation in (|52|) holds, then so does the second. The first relation, in terms of E(r), reduces 
to E' /E > (— 1/r) which is necessarily true. Hence we have obtained a class of solutions which show late 
time acceleration of the effective scale factor a-p, with the caveat that we must exclude a sphere around 
the origin from our domain of integration, the radius of this sphere being determined by the form of E{r). 

5 Discussion 

We have shown that it is possible within the framework of classical General Relativity, to construct models 
of universes in which the average behaviour of spatial slices is that of accelerating expansion. Although the 
LTB models are unrealistic (since they place us at the center of the Universe), they are useful in building 
intuition. Especially, since LTB is an exact solution, it helps towards a deeper understanding of averaged 
inhomogeneous cosmological solutions of Einstein equations. In particular, our solution could be assumed 
to apply not necessarily to the whole Universe, but only to a local underdense region such as a void. The 
volume average of the late Universe is dominated by voids, and structures occupy a tiny fraction of the 
volume. The average over such a distribution does not lead to a FRW model. The curvature of voids 
can be estimated to be proportional to minus the square of their Hubble expansion rate and thus must be 
negative [22]. The negative curvature LTB model discussed in this paper could be useful for describing 
such a situation. 

In both the examples which we gave, namely the power law models and the models given by Eqn. (|40|) . 
the qualitative behaviour of the evolution of the effective scale factor was independent of the size of the 
averaging domain. Whereas the power law models were pathological in the sense that, among other things, 
their spatial Ricci scalars diverged at the symmetry center, the second class of models described by (|4()jl 
had no such pathology. 

The point to be emphasized, though, is that our analysis clearly shows the importance of curvature 
(expressed as a non-zero energy function E{r)) in causing the late time acceleration. The fact that the 
limit of vanishing matter density (e — > 0) is well defined, means that the qualitative behaviour of curvature 
dominated models is not expected to change by adding a finite but small amount of matter. 

A few remarks comparing the results of the present paper with the concordance model in standard 
cosmology are in order. One could assert that observational data show that our currently accelerating 
universe has zero spatial curvature. On the other hand, our LTB model with zero spatial curvature (the 
marginally bound case) shows no acceleration. It might hence appear that our curvature dominated LTB 
model is of very limited interest. However, observations in the late Universe must not be matched only 
with a zero-curvature Universe. This assumption may be good in the early stages of evolution, but it 
is just a fitting ansatz to the late-time inhomogeneous Universe. What the consideration of averaged 
inhomogeneities in the present and other similar papers demonstrates is that the real Universe could in 
principle have negative curvature and yet exhibit acceleration. If this were indeed to be the case, it could 
eliminate the need for a dark energy. 

It might also be said that observational data show that the Universe has up to thirty percent of its 
content in luminous and dark matter and hence our low density curvature dominated model showing 
acceleration does not meet this criterion. However, in reality our model is in disagreement with the 
concordance model, which while being one of the simplest and most favoured possibilities, need not turn 
out to be the final answer. The consideration of averaged inhomogeneities shows that low density, curvature 
dominated models can also in principle fit the observational data. This issue should hence be regarded 
as an open one. Furthermore, as noted above, our LTB model could by itself be considered relevant for 
describing a locally underdense region such as a void. 

Our results highlight the intimate connection between the averaged spatial curvature, and the evolution 
of the kinematic back-reaction, as anticipated from the integrability condition given by Eqn. (|12|) . It 
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will be interesting to consider the more general dust models described by the metric and enquire if 
acceleration can again be produced in the approximation where the averaged three-curvature dominates 
over the averaged matter density during some epoch of cosmic evolution. 
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